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Abstract 

The interaction of the witness beam with the surrounding 
plasma particles and wakefields was studied. The impli¬ 
cations of the elastic scattering process on beam emittance 
and, emittance evolution under the focusing and acceleration 
provided by plasma wakefields were discussed. Simulations 
results from GEANT4 are presented in this paper. 

INTRODUCTION 

The next generation of particle physics colliders will need 
to supplement pp collisions with e + e~ and ep collisions 
to deliver precision and to address QCD research needs. 
Generally each successor collider should push the limits of 
the energy frontier further. 

Plasma accelerators have made tremendous progress in 
the last few decades since the inception of the idea from 
Tajima and Dawson j 1 ]. Nowadays, laser wakefield accelera¬ 
tors (LWFA) can achieve MeV-GeV level, electrons through 
millimetre to centimetre plasma cells ||2jQ. Electron beam 
driven plasma wakefield acceleration (PWFA) has demon¬ 
strated energy doubling for an ultra relativistic 42 GeV elec¬ 
tron beam in a metre long plasma structure [7). The acceler¬ 
ating gradients measured in these experiments can be in the 
range of 10-100 GeV/m, which are 3-4 orders of magnitude 
larger than that in today’s conventional RF-based particle 
accelerators. 

Towards the realisation of a collider scheme based on 
plasma wakefield acceleration, challenges and issues must 
be explored. 


EMITTANCE GROWTH DUE TO 
COULOMB SCATTERING 

Under the conditions where a beam travels in the vacuum 
with a constant acceleration, emittance decreases with in¬ 
creasing energy according to the conservation of the area 
in the phase space given by Louville’s theorem. This phe¬ 
nomenon is known as adiabatic damping. However, if the 
particles in the beam encounter a medium of gas or plasma, 
emittance diffusion occurs through scattering and competes 
against the adiabatic damping as suggested in EqJT] (8j ; 
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where 7/ is the rate of change of the acceleration, 7 f and 
7 i are the final and initial beam energies, respectively. F is 
written as, 
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where r e is the classical electron radius, m is the mass of 
the electron; n is the number of scattering centres, a 0 is 
the initial beam size interacting with the scattering medium. 
The constant accelerating field is given as E Z Qsin{(j)). 

Minimum and maximum scattering angles are determined 
through uncertainty in the momentum of the incident parti¬ 
cle, p, (Eq[3]) and the impact parameter, b. 


the quantum mechanical limit Ap b > h applies resulting in 
Eq|4j where h is the reduced Planck constant, 


J mm, max 


ph 


max,mm 


( 4 ) 


The maximum impact parameter, b max , comes from the 
shielding effect of the atomic electrons for a linear plasma 
wakefield. In a fully ionised plasma this will correspond to 
the Debye length, shown in Eq|5j where eo is the electric 
permitivity of vacuum, ks is the Boltzmann constant, T e is 
the temperature of the plasma electrons, n e is the number 
of electrons in the plasma and e is the charge of an electron: 
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The Debye length is the distance over which the potential of 
the nucleus is reduced to 1/e of its maximum value within a 
plasma due to the screening effect of the surrounding plasma 
electrons. For the nonlinear (bubble or blow-out) regime, 
where number density of the drive bunch is larger than the 
plasma, the maximum impact factor corresponds to the ra¬ 
dius of the ion cavity 

The minimum impact parameter can be related to the 
effective Coulomb radius of the nucleus, R. The extrema of 
the scattering angle can be rewritten as in Eq|6] 
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MONTE CARLO SIMULATIONS 


The above theory can be examined by comparing it with 
the results from a Monte Carlo code which can simulate the 












particle-matter interactions such as GEANT4 [ [T0|[TT| . A 
particular scenario was simulated where an electron beam 
with given parameters, under constant acceleration and fo¬ 
cusing, travels through a defined gas column undergoing 
only elastic Coulomb scattering. An example wake field of 
a 250 MeV proton drive beam is shown in Fig [T| within the 
simulation window of a few plasma periods. This result from 
LCODE [ 12) simulations suggests a 0.5 GV/m longitudinal 
field accompanied by a 0.1 GV/m radial field both of which 
were implemented in the GEANT4 model. 



Figure 1: Longitudinal and radial plasma wake fields due to 
a proton bunch that is located at zero. 


The initial beam consisting of 10 k electrons at 10 GeV 
was generated with randomly assigned positions and angles 
within Gaussian distributions with one standard deviation 
of 10 /im, 10 /irad, respectively, as shown in Fig. [2] These 
initial values were chosen considering a realistic emittance 
of an electron beam of 10 GeV. 

Particles were tracked 500 m through neutral Lithium (Li) 
gas (Z = 3, a = 6.941 g/mole) with a density of 6 x 
10 14 cm -3 . Li gas was chosen due to its orders of magnitude 
low scattering cross section compared to the other candidate 
media such as Rb (Z = 37). In reality, plasma is produced 
by ionisation of a channel through a chamber filled with a 
given gas with a radius given by the ionisation laser spec¬ 
ifications & Therefore particles travelling through the 
centre of the chamber may interact with the plasma ions and 
electrons as well as the surrounding neutral gas when they 
are scattered out of the plasma channel. With the current 
technology a plasma column on the order of a mm is possi¬ 
ble to produce. Nevertheless, in these simulation studies an 
arbitrary plasma radius of 100 mm was chosen to provide 
enough clearance for the particles. 

This 500 m plasma column with 100 mm radius was split 
into logical sections of 20 m, and the transverse phase space 
of the beam was reconstructed at the end of each section. 
For the analysis, primary particle tracks were isolated from 
possible secondary particles. 



Figure 2: Initial particle position (left) and angle (right) 
distributions provided for the simulations. 

Reconstruction of particle angles and phase space 

In order to produce a transverse phase space distribution 
after each section, angles of each primary particle should be 
determined at the end of the section under study. These can 
be calculated as mean angles corresponding to the last n scat¬ 
tering events as depicted in Fig. [3] Two example transverse 



Figure 3: The mean angle of each particle can be determined 
for the last n scattering events before the end of each section. 

phase space distributions, that are constructed from these 
tracking data, are presented in Figure [4] A geometric cut 
was implemented to reject the particle tracks falling outside 
of the plasma column with a radius of 100 mm. The rms 
emittance values were calculated (Eq|7} and Courant-Snyder 
parameters were extracted. The subsequent terms of EqjT] 
are implemented as in Eq[8] [9} [10| respectively, 
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where N is the number of particles considered in the cal¬ 
culation, Xi and x\ are the coordinates in the phase space. 
The evolution of the normalised emittance growth in the 
gaseous media and theoretical prediction are presented in 
Fig. [5] Green curve is the growth predicted in [8) whereas 
red curve represents the same model with the beam parame¬ 
ters used in this study. Theoretical emittance growth values 
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Figure 4: Reconstructed phase spaces at (a) 100 m and (b) 
500 m during the beam passage through a Li gas column. 
Particles within the lcr of the Courant-Snyder ellipses are 
shown in red. 


were calculated for the a Debye length of 0.3 fim assuming 
a 1 eV of excess energy for the plasma electrons at the given 
particle density. The GEANT4 results are given in blue 
dots. Accordingly, at 480m the emittance growth due to col¬ 
lisions is 1.1 x 10 -5 of the total emittance of 1.3 mmmrad. 
This is only a fraction of the total beam emittance. Another 
cause of emittance growth in a plasma accelerator can be 
mismatching between radial plasma focusing and the beam 
beta function. 


MATCHING BEAM AND PLASMA 

Beam envelope can undergo betatron oscillations due 
to the focusing component of the wakefields, E r . Beam 
emittance can be minimised when plasma focusing term, 
K = eE r /rm e ^c 2 , compensates against beam divergence 
by satisfying K « e 2 /cf q (l4, 15 ]. This oscillations are 
currently under study in the GEANT4 model to assess the 
access growth in emittance that might have been caused 
by the unmatched beam and plasma conditions, apart from 
the scattering contribution. The preliminary prediction is 
an initial emittance of 1 mm mrad and a plasma density of 
1 x 10 18 m -3 should provide the matching between beam 
and plasma and hence the minimum emittance. Studies to¬ 



Figure 5: The evolution of emittance growth throughout a 
gas column of 500 m according to the theoretical predictions 
and GEANT4 results. 


wards a numerical demonstration of the phenomenon are 
ongoing. 

CONCLUSIONS AND OUTLOOK 

As an advanced accelerating technique plasma wakeheld 
acceleration has ever-increasing prospects. Therefore, the 
potential issues of the scheme must be assessed carefully. 
This study was initiated to seek out the impact of the interac¬ 
tion of a witness beam with the surrounding plasma formed 
to provide acceleration. The emittance growth induced in 
the beam via beam-gas elastic scattering was studied, nu¬ 
merically. It has been analytically calculated for the given 
parameters that the beam-plasma interaction can cause an 
emittance growth of 0.014 nm after 500 m travel within the 
plasma due to multiple Coulomb scattering. In order to 
analyse the total beam emittance, the sensitivity of our cur¬ 
rent GEANT4 model for the unmatched beam and plasma 
parameters are under study. 
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